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ABSTRACT 


By its very nature, an optical heterodyne receiver is both a receiver and 
an antenna. Certain fundamental antenna properties of heterodyne receivers are 
described which set theoretical limits on the receiver sensitivity for the 
detection of coherent point sources, scattered light, and thermal radiation. 

In order to approach these limiting sensitivities, the geometry of the optical 
antenna-heterodyne receiver configuration must be carefully tailored to the 
intended application. The geometric factors which affect system sensitivity 
include the local oscillator (LO) amplitude distribution, mismatches between the 
signal and LO phasefronts, central obscurations of the optical antenna, and 
nonuniform mixer quantum efficiencies. The current state of knowledge in this 
area, which rests heavily on modern concepts of partial coherence, is reviewed. 

Following a discussion of noise processes in the heterodyne receiver and 
the manner in which sensitivity is increased through time integration of the 
detected signal, we derive an expression for the mean square signal current 
obtained by mixing a coherent local oscillator with a partially coherent, quasi- 
monochromatic source. We then demonstrate the manner in which the IF signal 
calculation can be transferred to any convenient plane in the optical front end 
of the receiver. Using these techniques, we obtain a relatively simple equation 
for the coherently detected signal from an extended incoherent source and apply 
it to the heterodyne detection of an extended thermal source and to the back- 
scatter lidar problem where the antenna patterns of both the transmitter beam 
and heterodyne receiver must be taken into account. Finally, we consider the 
detection of a coherent source and, in particular, a distant point source such 
as a star or laser transmitter in a long range heterodyne communications system. 
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1. INTRODUCTION 


Heterodyne or coherent detection can be advantageous in a variety of 
applications. Heterodyne receivers have at least two features which are quali- 
tatively different from incoherent (or direct detection) receivers (ref. 1) . 
First of all, the receiving bandwidth is determined by the IF bandwidth which, 
in principle, can be varied at will to give very high spectral resolution. 
Secondly, information related to the phase of the radiation signal is retained 
in the IF output and the outputs of two or more receivers can be correlated to 
make coherence measurements comparable to the aperture synthesis techniques of 
radio astronomy. 

To achieve high spectral radiation with incoherent or direct detection 
systems, radiation filters or spectrometers must be utilized and the combination 
of very narrow bandwidth and high sensitivity (low loss) is usually difficult to 
realize. In general, a heterodyne receiver will be more sensitive than a direct 
detection receiver with an equivalent noise equivalent power (NEP) for spectral 
resolutions below a cutoff bandwidth which depends on the NEP and the infrared 
wavelength (refs. 1, 2) . The sub-Doppler spectral resolution of heterodyne 
receivers can be exploited to study the molecular constituents and kinematics 
of remote sources yielding specific information such as altitude profiles of 
absolute abundance of the species, vertical temperature profiles, and wind 
velocities (ref. 3) . In detecting extraterrestrial thermal sources, the infor- 
mation is gathered by passive heterodyne spectrometers whereas, in our own atmo- 
sphere or in planetary atmospheres visited by spacecraft, active backscatter 
lidars can be employed. In contrast to the above applications where the radia- 
tion signal is totally incoherent or only partially coherent, the signal from 
the laser transmitter in a heterodyne communication system (ref. 4) is coherent 
except as modified by atmospheric effects (ref. 5) . This article attempts to 
present a unified theory of heterodyne receivers which addresses the optical 
design considerations for all of these applications. 

A representative heterodyne receiver is illustrated in Figure 1. Signal 
radiation is collected by an optical antenna and focused, along with a local 
oscillator beam, onto a square-law frequency mixer operating at the radiation 
frequency. The latter beams have center frequencies Vg and V L and powers 
Pg and P L . The two frequencies mix to give an output spectrum centered at the 
intermediate frequency Vjp = Vg - where Vjp is much smaller than the 

infrared frequencies Vg and V L and typically on the order of a GHz or less. 
The resulting signal current is amplified by an IF amplifier of bandwidth Bjp 
and rectified by a nominally square-law detector to give a current output pro- 
portional to the power in the IF. This is usually input to a low-frequency 
filter or integrating circuit to further enhance the spectral resolution and/or 
sensitivity and is then recorded. 

Although the present article will address most factors influencing the per- 
formance of the receiver in Figure 1, it will emphasize the design of the optica 
front end of the receiver for a variety of applications and, in particular, the 
manner in which the optical antenna geometry and local oscillator distribution 
affect system sensitivity. In Section 2 of this paper, we review the noise 
processes relevant to the IF signal and discuss the system signal-to-noise in 
the IF in terms of an as yet undefined mean square signal current. Section 3 
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briefly outlines the sensitivity improvement achieved by time integration tech- 
niques. In Section 4, we address the calculation of the mean square signal 
current in the mixer plane for a general, partially coherent, quasi-monochromatic 
source and, in Section 5, demonstrate the manner in which the IF signal calcula- 
tion can be transferred to any convenient plane in the optical front end of the 
receiver. In Section 6, we apply the general result to the specific problem of 
coherently detecting an extended incoherent source. The results of that section 
are then applied to the heterodyne detection of an extended thermal source in 
Section 7 and to the backscatter lidar problem in Section 8 and some useful 
design guidelines are generated. In Section 9, we apply the results of Sec- 
tion 4 to the detection of a spatially coherent source such as a laser trans- 
mitter in a heterodyne communications system or a distant point source such as 
a star. 


2. THE SIGNAL-TO-NOISE RATIO OF A HETERODYNE RECEIVER 

The power signal-to-noise ratio of a heterodyne receiver is a measure of 
its sensitivity since setting the ratio equal to one permits calculation of the 
noise equivalent power (NEP) . It is given, in most cases of interest, by 
(ref. 1) 



power 


^' i m 2 y 

<i s 2 > + <i-T 2 > + <iB 2 > + <ij 2 > + Oa 2 > 


( 2 . 1 ) 


o 

We will leave the calculation of the mean square signal current <i M } to later 
sections and limit our present discussion to the various noise terms in the 
denominator of Equation (2.1). 

The local oscillator induced shot noise, or quantum noise, <^ig 2 ^> is 
often the dominant noise if h\) >> KT B where T B is the equivalent blackbody 
temperature of a thermal source lying inside the antenna pattern of the 
receiver. Shot noise is due to fluctuations in the rate of arrival of LO pho- 
tons. If the LO power is much greater than the signal power, the mean square 
shot noise current is given by 


<is 2 > _ 2 ^ eB IF i DC 


2|3e 2 B 


IF 


hv 


Jj ^ d ^D B Q^D) I l(?D> 


( 2 . 2 ) 


where i D( ~. is the DC current generated by the LO, e is the electronic charge, 
B IF i s intermediate frequency bandwidth, and hv is the photon energy. 

The integrand contains the detector quantum efficiency T)g and the LO intensity 
I B which are assumed to vary over the plane of the detector defined by the two- 
dimensional coordinate r^,. The parameter 6 equals 1 for photoemissive mixers 
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while, for photoconductors, it equals 2 due to fluctuations in the generation 
and recombination of charge carriers as described by Levinstein (ref. 6) . 

One can rewrite Equation (2.2) in the more familiar form 


<i s 2 > » 


2f|gge B if P l 
hv 


(2.3) 


if we define an average quantum efficiency n 0 by 


n Q = 


IL 


dr D I L^ r D^ 


ff D 


(2.4) 


dr D I L^ r D^ 


and P B is the local oscillator power incident on the detector. 


Radiation from a thermal source contained within the receiver field of view 
and the receiver bandwidth B IF will be coherently detected and subject to 
so-called "heterodyne amplification." In some experiments, such as in passive 
heterodyne spectrometry, this thermal source is the object of study, while in 
others it corresponds to unwanted background noise. We will show in later 
sections that it can be described by the equation 


<i T 2 > 


ZPg^e b if P l 
hV(e hv/KT - 1) 


(2.5) 


where is an overall efficiency which depends in part on the design of the 

optical front end. 

Fluctuations in background radiation, which spectrally is outside the 
receiver bandwidth but within the infrared response band of the mixer, will also 
produce noise currents, given by <^i B 2 )> in Equation (2.1) as will sources of 
radiation outside the antenna pattern of the receiver but inside the heterodyne 
receiving bandwidth. McLean and Putley (ref. 7) have derived expressions for 
this noise component which are complicated functions of wavelength, spectral 
interval, detector area and temperature, and field of view. The latter noise 
is not amplified by the heterodyne process, however, and can be rendered 
negligible by choosing a large enough local oscillator power and by spatially 
and spectrally filtering the input radiation. 

Two other important sources of noise are Johnson or thermal noise asso- 
ciated with the mixer and the IF amplifier. The mixer noise is given by 
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<ij 2 > = 


( 2 . 6 ) 


4 kt m b if 


% 


where T M and Rj^ are the mixer's (or mixer load resistor's) temperature and 
resistance as seen by the IF amplifier. For most cooled mixers, this would be 
negligible compared with the amplifier noise given by 



4 kt a b if 

MR* 


( 2 . 7 ) 


where and R ft are the amplifier's noise temperature and input resistance, 

and M is a factor less than unity which accounts for impedance mismatches 
between the mixer and amplifier. 

Clearly, other sources of noise exist. "Excess noise" is common in 
receivers which employ diode laser local oscillators and generally arises from 
multimode effects or other non-ideal behavior in the LO. Noise can also be 
introduced at the electrical contacts to the mixer element or by temperature 
fluctuations in the mixer. These sources are unique to specific systems and 
will not be considered further here. 

With sufficient LO power, most of the above noise sources can be made 
negligible relative to the quantum noise (ig 2 ^> and/or the background thermal 

noise contribution ^i T 2 ^>. If the mean square signal current is given by an 
expression of the form 

2 

= 2 | ^Q 8 HET (hvy P S P L (2 ’ 8 


where Pg is the received signal power and ^HET 4s an as undefined 

heterodyne receiver efficiency, then, under strong LO illumination, the signal- 
to-noise ratio tends to 


< 4 M 2 > 


n HET P S 


N/ 


power <i 2 > + <i B 2 > 


hVBi F |s + f| T [exp (hv/KT) 


- r} 


( 2 . 9 ) 


Setting the latter ratio equal to one and solving for p s/ B IF yields the noise 
equivalent power per unit bandwidth; i.e.. 


NEP (W/Hz) 


hv 

rt HET 


•j B + r) T [exp (hv/KT) 



( 2 . 10 ) 
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where r) HET anc ^ Ht k°th depend on the optical front end geometry. In the 
quantum noise limit (hV >> KT) , Equation (2.10) reduces to 

NEP (W/Hz) = (2.11) 

n HET 


whereas, in the thermal limit (hV << KT) , it becomes 


n T 

NEP (W/Hz ) = KT (2.12) 

^HET 


If we include mixer and amplifier Johnson noise, we can write for a general 
photoconductor 


NEP (W/Hz) 


2hV 

4 - 

^HET 


h T hv 

n HET[ eXp(hV/KT) 



K ^ T M + V 

G 


where G is the "conversion gain" defined by Arams et al. (ref. 8). 


(2.13) 


3. DETECTION AND TIME INTEGRATION 

If the power signal-to-noise ratio in the IF is less than unity, the signal 
can be detected by integrating the detector output over a sufficiently long 
period of time. The voltage signal-to-noise ratio at the filter output in 
Figure 1 is linearly related to the power S/N by the equation (ref. 1) 



(*) J-(^IL] 1/2 

/power \l~2 / 


(3.1) 


The latter equation assumes that the IF amplifier has a rectangular bandpass 
spectrum (double sideband) , the rectifying detector is an ideal square-law 
device, the final output filter has a noise bandwidth B Q much less than B IF 
and the power S/N is much less than unity. Smith (ref. 9) has considered 
the more general case where the IF amplifier is not strictly square-law and does 
not have a rectangular bandpass spectrum. He has also considered power S/N 
ratios much greater than unity. If the output filter is a single stage RC 
circuit such that B Q = t 0 /4 = RC/4 , Equation (3.1) becomes 
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4. COHERENT DETECTION OF A GENERAL QUASI-MONOCHROMATIC SOURCE 

2 

We turn now to the calculation of the mean square signal current <(i M y 
for a general quasi-monochroma tic source. This problem has been considered pre- 
viously by Rye (ref. 10) and McGuire (ref. 11) . With only minor modification, 
the derivation given here parallels that of McGuire. If we assume that the 
detected radiation lies within a frequency bandwidth Avg that is narrow with 
respect to the center frequency Vg, the real signal field at the mixer plane 
can be represented by an expression of the form 

E s (r D ,t) = \{~2 Re{E s (? D ,t) e 10Jst } (4.1) 


where = 2TTVg and the complex signal field envelope E s (r D ,t) at the point 

r D in the detector plane varies slowly in time relative to the exponential 

exp(iajgt). The time dependence of the envelope might reflect the modulated 
output of a transmitter laser in a heterodyne communications system, the ampli- 
tude and phase fluctuations inherent in the signal from an incoherent thermal 
source or backscatter lidar, or even the effects of atmospheric turbulence on 
the signal. The envelope, through its dependence on the detector coordinate 
r D , also contains spatially dependent amplitude and phasefront information. 

If we represent the LO field by a similar expression, the current out of 
the square-law mixer is given by 


i (t) 




Eg (r D ) 




t) 


iw c t 
e b + 


^L ( r D ' 


ia) L t 


(4.2) 


where co L is the LO center frequency and the integral is over the active 
detector area. Upon performing the quadratic multiplication of fields in 
Equation (4.2), we obtain both sum and difference frequencies. High-frequency 
sum terms varying as exp (± 2 ic 0 gt) , exp (±2iu) L t) , exp (±i (ojg+co-^) t) , lie outside the 
bandwidth of the mixer and hence can be ignored. The difference terms produce 
two "DC" currents corresponding to the average signal and local oscillator 
induced currents and an additional mixing term given by 




(2e\ 

rr 

Uwj 

'J D 


^ * -+ iw IF t 

dr D r|g(r D ) Re<Eg(r D ,t) E L (r D ,t) e 


(4.3) 
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where the IF frequency tO IF = a)g - oj^. Squaring Equation (4.3) yields 


3-M 2 (t) 


2e 

hV 


dr D n Q ^ r D ^ 


x Re 


fl d?D ff B 

|e s ( r D ,t) E L *(r D ,t) e IF | Re|E s (r D ',t) E L *(r D ',t) e IF | 

2 (^) ff D d?D ff D d?D ' nQ(?D> 

x j^Re |E S (? D ,t) Eg* (? D 1 , t) E L (r D ',t) E L *(r D ,t)| 

+ Re | Eg (? D ,t) E s (r D ',t) E L *(r D ,t) E L *(r D ',t) e l2WlFt }J 


(4.4) 


If we average the above expression over a time interval T short compared to 
the coherence times of the signal and local oscillator field (T g and T L ) but 
long compared to the IF beat period, T IF , we may write 


2 i r t+T/2 2 

V ™ = T / dt (t) 

- , t-T/2 

~ ffv d ^ D ff n 2 (?D) 

X E g (r D ,t) Eg* (?d ’ , t) E L (r D ',t) E L *(r D ,t) (4.5) 


since the field envelopes can be viewed as effectively constant over this time 
interval and hence the terms varying as exp (±2iu)j F t) in Equation (4.4) average 
to zero over an IF beat period. In certain applications, such as passive 
heterodyne spectrometry of a thermal source, the integration time can be 
arbitrarily long. The limit of Equation (4.5) as T approaches infinity is 
then 
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/• t+T/2 

<i M > = lim - I dt i M 2 (t) 

T-*» ♦'t-T/Z 

= 2 (^) ff D d?D ^/' d?D ' n Q (?o) VV } 

x<E s (r D ,t) E s * (r D ' ,t)> <E L (r D ' ,t) E L *(r D ,t)> (4.6) 


where we have invoked the fact that the signal and local oscillator fields are 
statistically independent and hence the fourth-order correlation function 

<E S (r D , t) E s *(r D ',t) E L (r D ',t) E L *(r D ,t)^> can be written as the product of two 

second-order functions. The second-order correlation functions can be related 
to quantities appearing in the theory of partial coherence by noting that the 
"mutual coherence function" (MCF) of a quasi-monochromatic, stationary optical 
signal field is defined by (ref. 12) 


r 


s 



= <E s (r 1 ,t+x) E s *(r 2 ,t)> 


i w S T 


(4.7) 


Under the assumption of cross spectral purity (refs. 12, 13) , the spatial and 
time variables are separable leading to 


r s (r 1 ,r 2 ,T) = J s (^!'r 2 ) g ( t) e 1WsT (4.8) 

where g(0) = 1 and J s^ r l'^'2^ is t ' le " mutua l intensity function" (MIF) 
of the signal field. From Equations (4.7) and (4.8) we note that 
<Eg(r D ,t) E s *(r D ',t)>= Tg (r D ,r D ' ,0) = Jg(r D ,r D ’) and hence Equation (4.6) 
can be written in its final form 


2 

0-M y = 2 (hv) dr D J'J' D dr D n Q^ r D^ ^ r D ^ J S^ r D ,r D ^ J L^ r D ,r D^ 


(4.9) 


where and J L (r D ',rQ) are t ^ ie mu tual intensity functions of the 

signal and local oscillator fields in the detector plane. Calculation of the 
mean square mixing current by means of Equation (4.9) is not always a simple 
task due to the difficulty in computing Jg(r D ,r D ') for many sources of 
practical interest. In ensuing sections, we will demonstrate how the calcula- 
tion can be carried out in optical planes other than the detector plane and the 
enormous simplifications that often result. 
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Before closing this section, it is worthwhile to note two useful properties 
of the mutual intensity function; i.e.. 


J S* (r D' r D ,) = < E c*( r n' fc ) E C < r n' ' fc )> = J c ( r r,' ' r n) 


"S ' D S' D 


S D D' 


(4.10) 


and 


J S^ r D' r D^ I S 


(4.11) 


where ds t -'- me averaged signal intensity at the point r D . 


5. PROPAGATION OF THE MUTUAL INTENSITY FUNCTION 

Consider the signal electric field propagating from the antenna plane in 
Figure 2 to the detector plane. Small angle scalar diffraction theory (ref. 12) 
gives the electric field in the detector plane; i.e. , 


E s (? D ,t) 



P A (r A ) E s 


r A ,t 



ikr 


le 


AD 


Ar 


AD 


(5.1) 


where k = 2 it/A , P A<r A ) is the antenna pupil function and the term in brackets 
corresponds to a Huygen 1 s wavelet emanating from a point r A in the antenna 
plane and traveling a distance r^ to a point r D in the mixer plane. Then, 
from the definition of the mutual intensity function (MIF) , it is clear that 


J S ^ r D ' r D ' ^ -^E s (r D ,t) E s *(r D ',t)!> 






dr A' Pa^A 5 Pa^a') 


ik ( r AD _r A'D ' ^ 


^ r AD r A ' D ' 


E s[^A' t E S*(r a ' 't 


-A'D' 


"S \ A 


(5.2) 


For a stationary process, the time origin is of no consequence and therefore 


r 1+ . X AD r : i . 

E S l r A ' C j E S (a ft c 


r A’D’ 


E S ^ r A ' ^ E S 


r A " 


( r a 1 n 1 r in) 


(5.3) 
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Now, if the transverse dimensions of the antenna and detector pupil are small 
compared to the coherence length of the signal radiation defined by 1 = c/Av s , 


the variation of the signal electric field over a time interval 
t = (r A . D . - r p^) /c is negligible and Equation (5.3) is effectively the 
signal MIF in the antenna plane. Equation (5.2) then becomes the propagation 
law for the MIF as first derived by Zernike (refs. 12, 14); i.e.. 


J S ^ r D' r D ' ) 




P A (r A ) P A (r A ')| 


ik(r AD -r A 'D’ ) 


* r AD r A'D' 


J S ( r A' r A' > 


(5.4) 


If we substitute Equation (4.5) in (4.9) for the mean square signal current and 
reverse the order of integration, we obtain 


<iM 2 > 




dr A' p A< r A) p A ( r A ' ) J S ( r A ’ r A ' ^ 


X \^/y* dr D P D^ r D^ P D ( r D ' ^ J L ^ r D ' ' r D^ 


ik ( r AD "" r A ' D ' ) 


A r AD r A’D' 


(5.5) 


where P D (r D ) is the mixer pupil function. If we now define an effective 
local oscillator field given by 


^E^ r D ~ P Q ^ r D ' ^ E L^d' 

the corresponding effective MIF is then equal to 

J E^ r D ,,r D^ = J L^ r D ,,r D^ 


(5.6) 


(5.7) 


Substituting Equation (5.7) into (5.5) and comparing the resulting 
expression with the MIF propagation law (5.4), we note that the bracketed term 
in Equation (5.5) is simply the MIF of the effective local oscillator back- 
propagated to the antenna plane. We may therefore write for the mean square 
mixing current 


<:L m 2> = ffar A ffaZz 


P A< r A> P A< r A'> J S< r 


rj 1 ) Jp (rj 


A' L A 


,r A ) 


(5.8) 
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The physical significance of Equation (5.8) is that the calculation of 
mean square IF signal current can be carried out in any convenient optical plane 
as first pointed out by Rye (ref. 10) . This has practical importance since it 
is usually easier, for example, to compute the backpropagation of a coherent LO 
electric field through an optical system than to propagate the MIF of an 
incoherent source in a forward direction through the system to the mixer. This 
fact will be well illustrated in later sections. 

Although we have considered only free space propagation in the present 
derivation, the approach is equally valid when intervening optical elements 
such as lenses, mirrors, and apertures are present. The simple Huygens wavelet 
in Equation (5.1) is then replaced by an appropriate transmission function for 
the optical system (refs. 10, 12). 


6. HETERODYNE DETECTION OF AN EXTENDED INCOHERENT SOURCE 

The expressions derived up to this point have assumed a general, partially 
coherent, quasi-monochromatic source. We consider now an important practical 
application in which the signal radiation emanates from an extended incoherent 
source and propagates to the antenna plane as in Figure 3. The propagation of 
the MIF proceeds in precisely the same fashion as in the previous section 
except that there is no coherence between the Huygens wavelets emanating from 
the infinitesimal sources located at r s and r g ' . Thus the second-order 
correlation function in the source-antenna plane version of Equation (5.3) 
becomes 



(? s , t) 




<E s (r s ,t) E S *(r s ' ,t) > = I s (r s ) 


«(? s " 



( 6 . 1 ) 

where Ig(r s ) is the time averaged radiation intensity at the point rg in the 
source plane and 6(r g - r g ' ) is the two-dimensional Dirac delta function. It 
can be shown that substitution of Equation (6.1) into the source-antenna plane 
version of Equation (5.2) and performing the double integral over r g ' yields 
the propagation law for the MIF of an incoherent source (ref. 13); i.e.. 


J 


S 




I s (r s } 


, ik ( r SA- r SA') 

~2 

A r SA r SA' 


( 6 . 2 ) 


where the integral is over the finite dimensions of the source. We may now 
substitute Equation (6.2) into (5.8) and reverse the order of integration to 
obtain for the mean square IF signal current 
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dr s I s^ r s^ 



(6.3) 


Through use of the MIF propagation law given by Equation (5.4), we recognize 
the bracketed term in Equation (6.3) as the mutual intensity function of the 
backpropagated effective local oscillator (BPELO) evaluated at the points 
rg = r s ' . But, since J E (rg,rg) = Ip(rg), the time averaged intensity of the 
BPELO in the source plane, Equation (6.3) reduces to the relatively simple 
expression 



I S^ r S^ I E^ r S^ 


(6.4) 


Thus we have the very useful result that the mean square IF signal current is 
proportional to the overlap integral of the extended incoherent source intensity 
with the backpropagated effective LO intensity. In the next two sections, we 
will apply this result to the detection of thermal radiation and to the back- 
scatter lidar problem. 


7. THERMAL SOURCE DETECTION 

The total power AP radiated into a hemisphere, within the IF bandwidth 
Bj F , from a small area AA on a blackbody is 


at, 2tt 
AP = — 

A 2 


hVB 


IF 


j^exp (hV/KT) - lj 


AA 


(7.1) 


Only the power emitted in the direction of the receiver contributes to the 
signal MIF in the antenna plane. Thus, if the receiver is in a direction normal 
to the plane of the blackbody, we must multiply the above expression by a 
factor 1 /tt corresponding to the power emitted per steradian in the normal 
direction. We must also multiply by 1/2 to account for the fact that the 
heterodyne receiver detects only one polarization component. Thus, the 
intensity to be substituted into Equation (6.4) is given by 


_ II \/ l\ AP = hVB IF 

s = U AW AA = X 2 [ exp(hv/KT) _ f] 


(7.2) 
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and Equation (6.4) becomes 


Ze Dr-p /* /" -y -y 

<i M 2 >= p =r // dr s I E (r s ) (7 

hv[ex P (hv/KT) - 1] yy s b b b 

where the integral is simply the total backpropagated effective LO power sub- 
tended by the source. 

If the dominant noise mechanism is the LO-induced shot noise given by 
Equation (2.3), the IF signal-to-noise ratio is 


/S\ = < i i M 2 y 

) ^ /power y 6 j^exp (hv/KT) - lj 


where r| T is the overall heterodyne receiver efficiency for thermal source 
detection introduced in Equation (2.5) and defined by 


n T 



X E (r S> 


(7.5) 


where Lg is the average mixer quantum efficiency defined by Equation (2.4) 
and P L is the LO power incident on the detector. If the mixer quantum effi- 
ciency is uniform. Equation (7.5) reduces to 


Lip 




It. ( r c) 


(7.6) 


where we have used Equations (5.7) and (4.11). The quantity lL^ r S^ I s the 
intensity of the actual backpropagated LO rather than the effective LO. The 
quantity q T replaces the mixer efficiency in the corresponding equations in 
the classic paper by Siegman (ref. 15) . 


If the source is so large that the backpropagated LO is contained entirely 
within its disk radius, the integral in Equation (7.6) is simply the total LO 
power in the source plane. Except for an atmospheric transmission factor ri A , 
the latter is equal to the backpropagated LO power exiting from the antenna. 
Thus, the overall heterodyne efficiency (7.6) can be broken down into several 
components ; i . e . , 


n T 


VA n o n R 


(7.7) 


where q o takes into account routine optical losses due to reflections and 
scattering while ri R is a geometric efficiency which takes into account 
vignetting, central obstructions, LO phasefront curvature, etc. in the optical 
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antenna. Numerically, T| R is equal to the fraction of the original LO power 
which exits from the antenna during backpropagation . 

As an illustration, consider the Cassegrain telescope in Figure 4. Let us 
assume that the mixer is illuminated by the fundamental gaussian mode of the 
local oscillator laser. If the gaussian mode is not truncated too badly by the 
mixer perimeter or by the secondary mirror, it remains gaussian until it is 
truncated by the primary mirror of radius a and centrally obscured by the 
secondary of radius b in the antenna plane. The geometric efficiency r| R is 
then given by 


— 

/• a 2 

2 2 

2 


2 "J 

f dr re- 2(rAj) 
b 

= a 

- e~ a 

(7.8) 


where to is the gaussian spot radius in the antenna plane and we have defined 
two parameters (ref. 16) a = a/oo and y = b/a. The geometric efficiency has 
been plotted as a function of a and y in Figure 5. 

The important thing to note in Figure 5 is that, for a given nonzero value 
of the linear obscuration ratio y = b/a, the optimum efficiency is less than 
what one would expect based on simple blockage of the incoming radiation by the 
central obscuration. For example, y = 0.5 would imply an areal obscuration 
efficiency of 1 - y 2 or 75%. The peak efficiency in Figure 5, however, would 
only be about 47% if one were to choose an optimum gaussian spot radius corre- 
sponding to a = 1.3. Nonoptimum choices clearly result in significantly worse 
performance . 

Clearly, to maximize the efficiency of coherent detection of a thermal 
source which fills the receiver field of view, one wishes to choose an optical 
geometry which allows the effective backpropagated LO to exit from the telescope 
with near-unity efficiency. Although this is most easily accomplished with 
off-axis reflective telescope geometries which eliminate the central obscuration 
problem, one is not limited to such geometries in general. For example, if we 
use appropriate masks in the LO beam to create a local oscillator distribution 
in the mixer plane which matches the Airy pattern of the centrally obscured 
Cassegrain telescope in Figure 4, the backpropagated LO will form an annular 
disk in the antenna plane which matches the antenna pupil function and provides 
unity transmission. This result assumes, of course, that the mixer quantum 
efficiency is reasonably uniform. The transmission loss of the beam splitter 
in Figure 4 is included in the optical efficiency n o - 

For such large sources, the efficiency is not sensitive to the wavefront 
curvature of the LO beam except to the extent that it modifies the LO trans- 
mission through the antenna pupil. For example, if one considers two systems, 
projecting the same gaussian spot size in the antenna plane of Figure 4 but 
having two different radii of curvature for the LO phasefronts, the fractional 
transmission and hence the receiver efficiency will be the same. The system 
with the wider backpropagated LO divergence will detect point sources near the 
optic axis with less sensitivity but this will be compensated for by the 
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detection of additional point sources which are beyond the field of view of the 
receiver with the smaller backpropagated LO divergence. On the other hand, if 
the source is of limited spatial extent, maximum detection efficiency dictates 
that the backpropagated LO be contained totally within the source pupil function 
and hence LO phasefront curvature effects will play a more important role. For 
small thermal sources in the near field of the receiver, as in a laboratory 
experiment, this can be accomplished by choosing an optical system which effec- 
tively focuses the backpropagated LO onto the target source and provides near- 
unity transmission efficiency for the backpropagated LO. 


8. INCOHERENT BACKSCATTER LIDAR 

Consider the lidar system in Figure 6. An outgoing pulse of temporal 
width 6 is transmitted through the atmosphere illuminating the aerosol 
scatterers in its path. The mixer current at time t is due to radiation 
scattered at a time t - R/c from a volume defined by the length c6/2 within 
the receiver field of view as determined by the backpropagated effective LO 
intensity. Although the aerosol scatterers are randomly spaced and typically 
many wavelengths apart, the return is not strictly incoherent since the 
scatterers within the volume of interest are "frozen" in their positions during 
the passage of a short laser pulse, thereby producing a coherent or "speckle" 
component in the return. Thus, based on a single return, one cannot perform 
the long time average necessary to progress from Equation (4.5) to (4.6) in our 
derivation of the mean square mixing current However, if we imagine 

repeating the lidar experiment many times over the same source volume and 
obtaining an average current waveform out of the mixer, the coherent component 
would be expected to average to zero over the ensemble of measurements due to 
the random relative motions of the scatterers. After averaging a sufficiently 
large number of current waveforms, we would then be left with the incoherent 
component. Thus, if the physical process being observed is ergodic, i.e., 
ensemble averages are equal to time averages, the mean square mixing current 
will be given by where the notation now applies to either an ensemble 

average or time average since the two are equivalent. 

With the additional argument given above, we can apply Equation (6.4) to 
the pulsed backscatter lidar problem. The source intensity function Ig which 
is now a function of range (Z coordinate) as well as the transverse coordinates, 
is given by 

I S (R '^'S) s P p(R '^S>(ir) d ^ 1T ' ) I T ( R >r S ) (8-1) 

where I T (R,rg) is the intensity of the coherent transmitter beam at the 
range R and transverse coordinate rg , dCJ(Tr)/df2 is the differential scatter- 
ing cross section in the backward direction, c6/ 2 is the length of the 
scattering volume, p(R,r s ) is the density distribution of scatterers, and 
p is a factor of order unity or less which takes into account depolarization 
effects. The product [da (tt ) /dfT] I T (R,r s ) is the power scattered in the back- 
ward direction per steradian by a single scatterer located at the coordinates 
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(R,rg) while the product p(R,r s ) (cS/2) is the number of scatterers per unit 
cross-sectional area in the source volume. Substituting Equation (8.1) into 
(6.4) gives 


<iM 2 > = 



da (tt) 
dft 



P (R,r s ) Irp (R,rg) Xg(R,rg) 


( 8 . 2 ) 


which yields the important result that the mean square signal current is pro- 
portional to the overlap integral of three quantities - the coherent transmitter 
intensity, the backpropagated effective LO intensity, and the density distribu- 
tion of scatterers. It is useful to note that we have not made the assumption 
that the transmitted and local oscillator beams are coaxial in deriving Equa- 
tion (8.2). In fact, the equation can be used for bistatic lidar systems pro- 
vided the transmitter and receiver optical axes are nearly parallel and an 
appropriate offset between transmitter and LO beams is included before computing 
the integral. If the transverse separation between transmitter and receiver is 
small relative to the spot sizes of the transmitter and BPELO at the range R, 
the bistatic system can be treated as coaxial to a good approximation. 


As a simple numerical example, we now consider the case of gaussian trans- 
mitter and local oscillator beams described by 


and 



I L (R, r s ) 


2f l ‘ 2 U,(R>) 

2 e 
7TU) t (R) 

J_i 


(8.3) 


(8.4) 


where P T and P L are the transmitter and local oscillator output powers and 
w T (R) and (jl) l (R) are the corresponding guassian radii at the range R. Sub- 
stitution of Equations (8 . 3) and (8.4) into (8.2) yields 



2TT P [p(R)c5/2] d ^ - - ) - A 2 P t P l 
TT 2 CjO t 2 (R) + W l 2 (r[] 


(8.5) 


where we have assumed a uniform scattering density p (R) and a uniform mixer 
efficiency rig. Clearly, <^i M 2 ^ increases with decreasing W T and imply- 

ing that the signal level will be maximized in a laboratory scattering experi- 
ment by focusing the transmitter and backpropagated LO into the sample. 

If the scattering volume in the lidar system of Figure 6 lies in the far 
field of the transmitter and LO beam waists, we can use the approximations 
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oj t (R) ~ AR/trw TO and oi L (R) ~ Ar/ttoj^q where oj T q and u) L0 are the respective 
beam waists and R is the distance between the waists and the scattering 
volume (ref. 17). Equation (8.5) then becomes 


<i M 2 > = 2 


n Q e 

hV 


2wp[p (R) cS/2j 


da (it) 



( 8 . 6 ) 


• — 9 

which exhibits the familiar R dependence for the lidar equation. Equa- 
tions (8.3) and (8.4) suggest the definition of an effective area for the 

2 2 

gaussian beam waists given by A^, = /2 and A^ = /2. Further defin- 

ing an average antenna area A = (A,p + A^J/2 and letting A-^ = £A and 
Aip = (2 - £) A, Equation (8.6) becomes 





which has a maximum for £ = 1 given by 



(8.7) 


( 8 . 8 ) 


Thus, we have demonstrated that, if we constrain the sum of the transmitter and 
receiver areas to the value 2A, we obtain a maximum signal when £ = 1 or 
A l = Arp, i.e., when the antenna areas are matched. To include optical and 
atmospheric transmission losses. Equation (8.3) should be multiplied by 
and Equation (8.4) by E^Erq where rj A is the atmospheric transmission for 
the range R and P T q and TI^q are the efficiencies of the transmitter and 
receiver optical systems. 

It should be clear that, just as in the case of thermal source detection, 
any LO power falling on the mixer that cannot be backpropagated through the 
receiver optics to the source will contribute to the shot noise but not to the 
signal current and therefore represents a reduction in system signal-to-noise . 
Thus, vignetting, central obscurations, and phasefront errors can have a major 
impact on the lidar efficiency by (1) reducing the transmission of the back 
propagated LO and (2) influencing the antenna pattern of the backpropagated 
effective LO in Equation (8.2). The antenna patterns of vignetted, centrally 
obscured, and decollimated gaussian beams have been computed by Klein and 
Degnan (ref. 16) . 
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9. COHERENT SOURCE DETECTION 


For a spatially coherent source such as a laser or distant star, we can 
write for the mutual intensity function at the mixer 


J S ( r D ' r D ' ) 


[« 


(r D ) 



£ S ( r D’ ) 



(9.1) 


where £ s and <p s are real functions which describe the signal amplitude dis- 
tribution and phasefront in the mixer plane. A similar expression can be 
written for the laser LO. Substituting Equation (9.1) and the LO equivalent 
into our general expression for given by Equation (4.9), we obtain for 

a coherent source 




2 \ _ 


— V 

hv ) 


fl 


dr D 1 V r D ) 


e s (r D ) 


£ L (r D ) e 


i [jf>s ( r D^ “^L ( r D^ 


(9.2) 


In the trivial case where the mixer efficiency and the signal and LO beams are 
uniform over the mixer of area A D , Equation (9.2) reduces to the familiar form 


<iM 2 > = 


UQe 

hv 


P P 
S L 


(9.3) 


2 

where Pg = Eg A D . In the most general case, we can use Equation (2.8) to 
define a coherent heterodyne efficiency given by 


'HET 


n Q P L P S 


fl 


dr r 


,(r T 


>) £ s<^d) e l (£ d) ex P ai«f> s (r D ) 


■H* 



(9.4) 


where Pg is the total signal power in the mixer plane. Equation (9.4) can 
also be written in the form 


J 

fj d r D n g (r D ) e s (r D ) E L (r D ) expji 

[*S ^ r D^ “ ^L ^ r D^j | 

2 


JJ dr D r| Q(f D ) £ l (r D ) 

[fj 

f* 2 •*> 

dr D £ S 



where we have used Equation (2.4) and the explicit expression for P . Degnan 
and Klein (ref. 18) have performed computations of ^het ^ or case where the 

signal and LO phasefront curvatures are matched, the mixer efficiency fig is 
uniform and e s (r D ) is an Airy pattern formed by a centrally obscured, circular 
antenna illuminated by a plane wave from a distant point source. As the size of 
the central obscuration is increased, more incoming radiation is blocked by the 
obscuration and a smaller fraction of the radiation which reaches the mixer 
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plane is contained in the central lobe of the signal Airy pattern. Degnan and 
Klein (ref. 18) considered several illumination profiles for the LO including 
uniform, gaussian, and an Airy pattern matched to the signal Airy pattern. 

Their results are summarized in Figure 7. Optimum detection efficiency is 
achieved when the mixer captures the entire signal Airy pattern and a matched 
LO is used. In this instance, the receiver efficiency is simply 1 - y 2 
(where y is the obscuration ratio defined previously for the Cassegrain 
antenna in Figure 4) corresponding to the areal obscuration loss and repre- 
sented by the "matched" LO curve in Figure 7. The difference between the ideal 
or "matched" LO curve and the uniform or gaussian curves corresponds to the 
heterodyne detection efficiency ^het " 

If the mixer is illuminated by a uniform LO, the optimum Airy disk radius 
(to the first null) is found to be ~ 1.35R]-, where R D is the mixer radius. 

It should be noted that the Airy disk radius varies with the obscuration ratio 
for an optical antenna having a given f-number (ref. 18). The optimum effi- 
ciency Lhet as a PP rox; *- ma tely 83% for no obscuration and falls rapidly as the 
obscuration ratio is increased even if one chooses an optimum signal spot size. 
An optimized gaussian LO with waist radius 0 ) = 0.64R A and a central Airy 
signal disk which matches the mixer radius R^ yields greater sensitivity com- 
pared to the uniform LO since it more closely matches the intensity distribution 
of the central Airy disk for the signal. The power contained in the outer rings 
of the Airy pattern is lost, however, and this accounts for the major difference 
between the "ideal" matched LO and gaussian LO curves in Figure 7. For a more 
detailed discussion, and for more general plots of non-optimized geometries, 
the reader is referred to the original paper by Degnan and Klein (ref. 18) . 

It is a simple matter to compute the effects of misalignment between the 
signal and LO beams or of a mismatch between phasefront curvatures using the 
general expression (9.5). For example, if the two wavefronts are misaligned by 
an angle 9 in the y D direction illustrated in Figure 2, the exponential 
argument in Equation (9.5) is 


^S (r D } ' 4>L (r D> 


— (kc — 


kr) • r 


D 


k sin 


VD 


where kg and k L are the propagation vectors for the signal and LO beams, 
] kg | ~ |k L | ~ k = 2tt/A, and y D is the y-component of the vector r D - For 

cylindrically symmetric fields, Equation (9.5) reduces to a special case 
previously derived by Cohen (ref. 19); i.e.. 


fl HET 



r D 


£ S ^ r D^ 


e L (r D ) J Q (kr D sin 0) 



dr D r D 



( r D ) 



£ L 


(9.6) 
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where r Q is the radius of the mixer, and we have used y D = r D cos <|) D and 
the integral expression for the Bessel function J 0 (z), i.e.. 


J Q (z) 



iz 


e 


sin cj> 


(9.7) 


Cohen (ref. 19) has generated plots of r) HET for a variety of source-LO 
illumination function combinations such as uniform-uniform. Airy-uniform, 
matched Airy-Airy, uniform-gaussian, and Airy-gaussian . He considered the 
tolerance of the various combinations to misalignment and allowed for a quad- 
ratically varying mixer quantum efficiency. The sensitivity to misalignment 
for the various combinations varied less than 15% relative to the most sensitive 
uniform-uniform case given by 


n HET n Q 


2J Q (kr 0 sin 0) 
kr Q sin 0 


(9.8) 


Thus, 0HET = 0Q for no m isalignment and HhET = Bg/2 for 0 = 0.5A/ (2r 0 ) 
corresponding to a half-wavelength phase difference over the mixer diameter 2r 0 . 
For a wavelength of 10 ym and a mixer diameter of 200 ym, the misalignment angle 
at which the detection efficiency is reduced by a factor of 2 is 0 = 1.4°. 

For a mismatch in phasefront curvatures, the exponential argument in 
Equation (9.5) is 



4>r.(r D ) 



where C s and C E are the curvatures of the signal and LO phasefronts at the 
mixer plane. For cylindrically symmetric beams. Equation (9.5) reduces to 


^HET 



ng(r D ) e s^ r D^ 


e L (r D ) 
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r D e L 


V r D> 



dr D r D £ S (r D } 


(9.9) 
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For the uniform-uniform case 


sin 

^HET n Q \ 


and n HET = n Q for = (c^T " C^) = 0 While n HET = 0 for A (^) = 2X / r Q 2 

where r Q is the mixer radius. Thus, if the local oscillator beam has a 
planar phasefront (C E = 00 ) , the signal beam phasefront curvature must satisfy 

C S >> r o 2 / 2X - 

It should be noted in closing that we have arbitrarily chosen to perform 
the above calculations in the mixer plane. For a particular antenna or LO 
geometry, it may be more convenient to perform the computation in some other 
optical plane as noted previously in Section 5. 



10. CONCLUDING REMARKS 

This article has attempted to present a unified approach to the calculation 
of signal-to-noise ratios in optical heterodyne receivers for a variety of 
important applications. No attempt has been made to give an exhaustive review 
of the existing literature. The references cited are those which, in the 
author's opinion, either lend themselves particularly well to the development 
of the general theory of optical heterodyne receivers given here or have pre- 
sented numerical results having widespread application. There are, for example, 
various uncited articles which present calculations of signal-to-noise for very 
specific incoherent source or backscatter lidar geometries. These have usually 
employed brute force computational methods that give little insight into the 
general approach for optimizing system sensitivity. While these provide 
excellent tests of the general theory, the articles were deemed to be too 
specialized to be included in the present review. 

Clearly, no attention has been paid to the effects of the atmosphere on 
coherent wave propagation. Although the amplitude and phase fluctuations pro- 
duced by the atmosphere are inherently included in the complex electric field 
envelopes introduced in Section 4, no attempt has been made here to give a 
quantitative assessment of their impact. In the approach taken here, the atmo- 
sphere can be viewed as simply another optical element through which the coher- 
ent backpropagated effective LO must pass to reach the signal source or vice 
versa. In the thermal source detection and backscatter lidar problem, the atmo- 
sphere presumably modifies the backpropagated effective LO intensity distribu- 
tion thereby influencing the overlap integral in Equation (6.4). A number of 
papers in this area have appeared since the early work of Fried (ref. 5) 
including a rather extensive recent report by Capron et al. (ref. 20) appli- 
cable to coherent optical radar. 


482 



REFERENCES 


1. T.G. Blaney, Space Science Reviews ,17, 691 (1975). 

2. J.H. McElroy, Applied Optics, 11,1619(1 972). 

3. M.J. Mumma, T. Kostiuk, and D. Buhl, Optical Engineering, 
V7,50( 1978) . 

4. J.H. McElroy, N. McAvoy, E.H. Johnson, J.J. Degnan, F.E. 
Goodwin, D.M. Henderson, T. A. Nussmeier, L.S. Stokes. 

B.J. Peyton, and T. Flattau, Proc. IEEE, 65, 221 ( 1977) . 

5. D.L. Fried, Proc. IEEE, 55,57( 1967). 

6. H. Levinstein, Applied Optics , 4 , 639 ( 1965) . 

7. T.P. McLean and E.H. Putley, RRE Journal , 52, 5( 1965) . 

8. F.R. Arams, E.W. Sard, B.J. Peyton, and F.P. Pace, IEEE JQE, 
Qe-3, 11(1967). 

9. R.A. Smith, Proc. IEEE,98,43( 1951 ) . 

10. B.J. Rye, Applied Optics, 18, 1390( 1979) . 

11. D. McGuire, Optics Letters ,ji, 73 ( 1980 ) . 

12. M.Born and E.Wolf, "Principles of Optics", 5th Ed., 

Chapt. 10 (Pergamon,New York, 1975) . 

13. L.Mandel and E.Wolf, Rev. Mod. Phys . ,37,231 (1965) . 

14. F.Zernike,Physika,_5,791 ( 1938) . 

15. A. E. Siegman, Applied Optics ,j>, 1 588( 1966) . 

16. B.J. Klein and J.J. Degnan, Applied Optics , 13,2134( 1974) . 

17. A. E. Siegman, "An Introduction to Lasers and Masers", 

Chapter 8 ( McGraw-Hill , New York, 1971). 

18. J.J. Degnan and B.J. Klein, Applied Optics, 13,2397( 1974) ; 
Erratum, _L3 , 2762 ( 1974) . 

19.S.C. Cohen, Applied Optics, 14, 1953( 1975) . 

20.B.A.Capron,R.C. Harney, and J.H. Shapiro, "Turbulence 
Effects on the Receiver Operating Characteristics of a 
Heterodyne Reception Optical Radar", Project Report 
TsT-33 , Lincoln Laboratory. Massachusetts Institute 
of Technology( 1979) . 



BANDWIDTH BANDWIDTH 



Figure 1.- Block diagram of a representative heterodyne receiver. 



Figure 2.- Huygen's wavelet model for propagation of the mutual 
intensity function between the antenna and mixer plane. 



Figure 3 .- Huygen's wavelet model for propagation of the mutual 
intensity from an extended incoherent source. 
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GEOMETRIC RECEIVER EFFICIENCY, n R 












Figure 6.- Functional diagram of a heterodyne incoherent backscatter 

lidar system. 



Figure 7.- Maximum receiver efficiency factors in dB for detection of 
a distant point source by a heterodyne receiver consisting of a 
general centrally obscured telescope (primary radius a, secondary 
radius b) as a function of linear obscuration ratio y = b/a and 
several optimized LO distributions (uniform, gaussian, and matched 
Airy) . 
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